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Abstract. We consider a class of linear differential operators acting on vector- valued function 
spaces with general coupled boundary conditions. Unlike in the more usual case of so-called 
quantum graphs, the boundary conditions can be nonlinear. After introducing a suitable Lya- 
punov function we prove well-posedness and invariance results for the corresponding nonlinear 
diffusion problem. 



1. Introduction 

Following earlier intuitions, K. Ruedenberg and C. Scherr developed in 1953 a new technique 
with the aim of studying electronic properties of conjugated bond systems, and in particular of 
aromatic molecules in [RS53] . Their idea was to set up a Schrodinger equation acting on a quasi- 
1-dimensional domain that can be schematised as a network of atoms, the network's edges being 
the chemical bonds. Their results aroused broad interest in the community of quantum chemists 
and marked the birth of the so-called free- electron network model. 

Over two decades later, when most solid state physicists had already turned back to the origi- 
nal, computationally more feasible discrete tight binding approximations, the free-electron model 
began to be studied by analysts and theoretical physicists. Among the first results we mention 
well-posedness results for the heat equations on networks obtained by G. Lumer in |Lum80| . Fur- 
ther interesting results followed soon: among others, F. Ali Mehmeti, J. von Below, P. Exner, 
S. Nicaise, Yu. V. Pokornyl and J. -P. Roth extended Lumer's results considering more and more 
general node conditions, providing interesting descriptions of the spectrum, discussing nonlin- 
ear and/or higher dimensional problems and establishing an interplay with quantum physics and 
theoretical mechanics. We refer to |PB04[ IKuc08j for a survey of these early investigations. 

This topic has finally gone mainstream in the late 1990s, when T. Kottos and U. Smilansky 
have observed in jKS97[ IKS99b| that models based on differential operators on metric graphs can 
play a fundamental role in the theory of quantum chaos. Ever since, network-based differential 
models for Schrodinger equations (and, by extension, also diffusion, Dirac and Pauli ones) have 
been commonly referred to as "quantum graphs" in the literature. 

Another interesting development has begun with KS99a, IHarOOl IKuc04[ IPan05j . where a va- 
riety of nonstandard boundary conditions for partial differential equations on graphs have been 
described. In this general sense, a quantum graph is just a particular way to look at a vector- 
valued diffusion equation with coupled boundary conditions, and the associated elliptic problems 
can be studied by means of classical Sturm-Liouville theory: observe that the above formalisms 
contained in the above papers are essentially just simple cases of the general framework intro- 
duced in SS65, SS66 . Though, such representations of coupled boundary conditions have paved 
the road for the development of abstract functional analytical methods for the treatment of gen- 
eral diffusion problems. It has been observed by several authors that quantum graphs represent a 
handy source of examples for unusual or even pathological behaviours of a diffusion equation. 

Nonlinear Robin-type boundary conditions for parabolic equations on domains are relatively 
common in the literature. They model an outgoing flow that depends nonlinearly on the temper- 
ature or the density at the domain's boundary. However, they have been seldom considered in 
the framework of quantum graphs. To the best of our knowledge, in the specific case of networks 
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they have been treated only in |Bel91[ IBel93j . The more general case of differential inclusions 
on ramified spaces is treated in jAN93j . while nonautonomous semilinear parabolic systems have 
been treated by several autors following the seminal article [5ma88 • 

The aim of this note is to introduce a general theory of 1-dimensional parabolic systems fea- 
turing nonlinear boundary conditions. Thereto we apply the theory of gradient systems and a 
formulation of (coupled) boundary conditions that is strongly inspired by the aforementioned 
works by P. Kuchment. In Section 2 we are going to prove well-posedness for a large class of 
quantum graph-like systems with nonlinear boundary conditions in variational form. In Section 3 
we are going to discuss some of their qualitative properties and L p -well-posedness. We conclude 
the paper by briefly discussing the issue of diffusion equations with nonlinear dynamic boundary 
conditions in Section 4. 



Finally, we emphasise that while the theory of quantum graphs has been the source of inspiration 
for the present investigation, it is misleading to connect our results to those in that field. In fact, 
while in the linear case well-posedness and qualitative properties for diffusion and Schrodinger 
equations on networks (as well as on domains) often come in pairs, in the nonlinear case most such 
connections fail to hold. The reason is that, of course, spectral methods lose their strength and have 
to be replaced by a hard analysis approach which is specific to the considered class of differential 
equations. In fact, our approach is based on the theory of gradient systems, which is currently 
pretty much bound to real Banach spaces. This seems to prevent any investigation of Schrodinger 
equations, even of damped ones. This is why we prefer to refer to previous investigations on 
diffusion equations with nonlinear coupled boundary conditions, although Schrodinger equations 
with concentrated nonlinearities that may look formally similar to ours have also been treated in 
the literature, see e.g. |AT01j . 

2. General setting and well-posedness results 

Let H be a separable, real Hilbert space and Y be a closed subspace of H x H . Let T > 0. 
Throughout this section we discuss the vector-valued diffusion equation 



(AV) 



§- t u{t,x) = £su(t,x) + i>(t,x), t G [0, T], x G (0, 1), 

«(*)|{o,i} e Y, t£[0,T], 

Ma + ^((u(% ,i } ) e Y\ te[0,T], 

u(0) = Uq, 



for an unknown u : [0,T] x (0, 1) — > H with inhomogenous term ip : [0,T] x (0,1) — > H and a 
possibly nonlinear function <fi :Y —> Y, occurring in the boundary condition. Here we denote 

f\{ 0< i}-[ m ) and 

for any / : [0, 1] — > H smooth enough. If <p is a linear mapping and H is finite dimensional, it has 
been shown in Kuc04 that the boundary conditions appearing in (AV) are the most general ones 
leading to self-adjoint diffusion operators. 

Example 2.1. One of easiest nontrivial parabolic systems that can be written in the form (AV) 
is possibly that considered in [ZZ03 . Restricting for the sake of simplicity to the 1-dimensional 
case, their boundary conditions for the unknowns u\ , U2 read 

u'S) = (-l) l + 1 Aie pitllW+ « 1U2W , u' 2 {i) = (-l) l + 1 A 2 e P2 " lW+ « 2U2(i) , i = 0, 1, 

for some positive parameters Xi,Pi,qi, i — 1,2. Clearly, it can be rewritten as 

letting H := R 2 , Y := {0} and <fi := ((f>i,(f>2), where 



zA AieP l2l+9lZ2 



z 2 ' V \Z2 ' \X 2 e P2Zl+q2Z2 
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A thorough blow-up analysis of this system has been performed in [ZZ03] and in many subsequent 
papers. 

Example 2.2. Let Y — Range/, where I is the n x m (signed) incidence matrix of a finite, 
simple, directed graph with node set V = {v 1; . . . , v„} and edge set E = {e 1; . . . , e m }, I + ,I~ are 
the matrices whose entries are the positive and negative parts of the entries of I and 

(2.1) /:= ' 1 > 



If <f> = 0, then the boundary conditions of (AV) agree with conditions of continuity /Kirchhoff- 
type, cf. [KMS07] . More generally, if <f> = ID^ 1 B(ID~ 1 ) T , they agree with so-called 5-coupling 
conditions 

u e (t,v) = u f (t,\/) =: u v (t), (e[0,T],e~v~f,veV, 

E e ~v^(M) = £ w6V *w«w(*), te[o,r|, vev, 

associated with an n x n-matrix B — (b vw ), cf. Mug07 . Here we write e ~ v if the edge e is 
incident in the node v, D is the diagonal matrix of vertice degrees and B is a n x n-matrix. 

Example 2.3. In the case of H = R m and 

Y := (1) x H = {(c, c, . . . , c,ai,a 2 , . ■ . ,a m ) e H x H : ai, a 2 , . . . , a m ,c G R}, 

(AV) reduces to a diffusion problem on a metric star with m edges with nonlinear S-coupling 
conditions in the origin and Neumann boundary conditions in the external nodes. 

For m — 2, the arising boundary conditions agree with those considered in [ATOlj . There, 
H 1 -well-posedness for a Schrddinger (instead of heat) equation and concentrated nonlinearity (at- 
tractive or weakly repulsive interactions ) in the origin has been proved. 

A problem similar to (AV) has been discussed also in [CGlOj . as the authors search for the 
viscosity solution of a hyperbolic system that arises in the Lighthill-Whitham traffic model. They 
show that a semilinear heat equation with concave nonlinearity on a star has a unique solution, 
for all small initial data in H . 

We incorporate the first boundary condition of (AV) in the Banach spaces 

H^:={feH x (0,l;H):f l{0;1} €Y} resp. Cy := {/ G C([0, 1];H) : /|{o,i} G ^} • 

For the inner product and the norm in H Y , inherited from H 1 ^, 1; H), we write (■ , ■)#!. resp. 
|| • \\ H i . For the norm in Cy, inherited from C([0, 1];H), we write || • \\c Y - We also write 

X 2 :=L 2 (0,1;H) and \\ ■ \\ 2 := \\ ■ \\l*(o,1;H)- 

Our aim is to re-write (AV) as a nonlinear abstract Cauchy problem - more precisely, as a 
gradient system. To begin with, we assume that <f> = V$ for a continuously differentiable function 
$ : Y — » K and introduce the functionals 

£o(f):=lj o \\f'(x)\\ 2 H dx, feH Y , 
Siif) ■= Hf\{o,i}), f G C Y , 

and 

(2.2) £ := £ + £i : H Y R. 

Lemma 2.4. The Hilbert space H Y is densely and continuously embedded into X 2 . 

Proof. The assertion follows from the inclusions Hq(0, 1; H) C H Y C i/ 1 (0, 1; H), since Hq(0, 1; H) 
is densely and H 1 (Q, 1; H) is continuously embedded into X%. □ 

Remark 2.5. We observe that 

(2-3) ||/| {0 ,i}lkxH < ^WfU < C\\f\\ H , , feH Y , 

for some C > 0. The former inequality follows from 

ll/|{0,i}l| 2 ff x ff = ll/(0)|| 2 ff + ||/(l)|| 2 ff <2||/|| 2 oo , 
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and the latter from the continuous embedding o/iJ 1 (0, 1; H) into C(0, 1;H). 

Inequality (|2.3p implies that £ : H Y — > K defined in (|2.2j) is a continuously diffcrcntiablc 
function. We recall that the gradient V£ with respect to the Hilbert spaces H Y and X 2 is the 
(generally nonlinear) operator V£ : D(V£) — > X 2 given by 

D(V£) :={/£^:3(6l 2 s.t. £'(f)9 = (£ , 9) x . 2 V9 e H Y ), 
V£(/) :=£, 

where £' : -ffy — > (H Y )' denotes the derivative of £. 

Next we are going to show that the initial-boundary value problem (AV) corresponds to the 
gradient system below with respect to H Y and X 2 

f u(t) + V£(u(t)) = *(t), te[0,T], 
\ w(0) = ito, 

where := ^>(i, •)• 

To this end, we observe first that the derivative of £ is given by 

£'(/)*= f 1 (f(x) , 6'(x)) H dx + (4>(f mi} ) , 9 mi} ) H xH, f,6€H Y . 
Jo 

Take / G D(V£) C and choose 6 e H%(0, l;H). Then, by definition of D{V£), there exists 
such that 

£'(f)0= (f(x),6'(x)) H dxl (Z(x),6(x)) H dx. 



We infer £ = — /" and / e -£f 2 (0, 1; H). For arbitrary 6* € _ffy we obtain, applying integration by 
parts, 

£'(f)0 = - J {f"i*) > + . #|{0,1}>™ + (0(/|{O,l}) . %o,i}>™ 

, dx = - [ (f'(x) , 0{x)) H dx. 

Jo 

Since one can find for every y £ Y a 9 € fly with #|{o,i} = V, h follows that ^ + </>(/|{o,i}) € 1^. 

If, on the other hand, u(t,-) € i?y n H 2 (0,1; H), t <E [0,T] satisfies the boundary condition 

a«(t) 



9k 



</> (w(t)|{ .i}) € F^, the above calculation yields € D(V£). 



The next Lemma is of rather general nature and similar to Mug08, Lemma 2.1]. We recall some 
definitions and well-known facts: Let V, X be Hilbert spaces such that V X, i.e. V is densely 
and continuously embedded into X. A function £ : V — > R is called coercive if for every c € K the 
sublevel set {/eF: £ (/) < c} is bounded in V. It is called X-elliptic if £"(•) := £ (•) • ||^ is 
convex and coercive for some w > 0. Obviously £ is coercive if £{■) > a\\ ■ \\ Y — (3 for some a > 
and /3 > 0. For quadratic forms the converse also holds true with (3 = 0. 

Lemma 2.6. Let V,A &e as above. Let £$ : V — >• M fee a quadratic form and £\ : V — > R &e 

convex and satisfy 

(2.4) > -fc||/||x e -j8 for all f eV and some (3,k>0, 

where X e is some Banach space such that V X e «-»• X and verifying the interpolation inequality 

||/lk <M £ ||/|| y ||/||^ £ , /GV, 

/or some e €E [0, 1) and some M e > 0. Then £q + £1 is X-elliptic if £q is X-elliptic. 
Condition (|2.4p is m particular satisfied if additionally £\ is continuous. 
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Proof. First we estimate E\ by 

£i(/)>-fe||/lli-/? 

>-^||/||y||/|U e -/3 

>-kKM e \\ff+*\\ftx e -f} 

> -kKM e (S\\f\\ 2 v + C(S)\\ff x )-p 

where K is a constant resulting from the continuous embedding of V into X e . The last estimation 
follows applying Young's inequality 

xy < Sx p + C(S)y p ', x,y>0, 

with p = 2/(1 + e) and 1/p+l/p' = 1. 

Assume now £o to be A-elliptic. Let / G V. Then for appropriate a > and w > we have 

£>(/) +u\\ffx >a\\f\\ 2 v . 
Choosing S small enough, we obtain 

(£o + &)(/) >a'||/|| 2 y-^||/||i-/3 

where < a' := a — kKM e 5 and < u/ := oj + kKM e C(5). Furthermore, (£q + is convex if 
£q is convex. 

The last assertion results from the Hahn-Banach theorem or directly as following. Let £i(0) = 
—6 and e > 0. Then there exists a 8 > 0, such that £i(f) > — b — e for all / with ||/|| < 8. Now 
convexity of E\ yields for every / with ||/|| > 6 

^cfl+a-^m^/)*-*-. 

and together with £i(0) = —b we get £ i(/) > — 1||/|| — b whenever ||/|| > S. Hence, we can choose 
7 = b + e and k = 4 to obtain 

> — ^||/||x E — 7 for all / € V and some 7, k > 0, 
which clearly implies (|2.4p . □ 

Now we can state the main result of this section. 

Theorem 2.7. Let (f> : Y — > Y" be a function that satisfies (f> — V$ for some convex $ € 
C (Y, R) and that maps bounded sets into bounded sets (what it does automatically ifY is finite 
dimensional). Then for all initial data uq £ Hy and all "J" € L 2 (0,T;^2) there exists a unique 
solution u e H 1 (0,T;X 2 ) n L°°(0,T;4) to (AV). 1/ * = 0, t/ierc tte solution also belongs 
to L 2 (0,T;H 2 {Q,l;H)) for every T > and even to L 2 (0, 00; H 2 (0, 1; H)) if in addition $ and 
therewith £ is bounded from below. 

These strong regularity properties of the solution are the main advantage of the approach 
based on gradient systems over the more general one based on maximal monotone operators as 
in |AN93[ IFGGR06] . For the proof of Theorem 12.71 we need the following general result based 
on [Llo69l Thm. 2.1.2 bis, p. 163], see also [CFlOl Thm. 8.10. 

Theorem 2.8. Suppose that V is a reflexive, separable real Banach space, which is densely and 
continuously embedded into the separable, real Hilbert space X. Suppose that £ : V — > R is a 
continously differ entiable function such that the derivative £' : V — > V maps bounded sets into 
bounded sets. Assume that £ is X-elliptic. Then, for all T € (0, 00), all \& G L 2 (0,T; X) and all 
initial data Uq € V there exists a unique solution u G H l (0, T; X) n L°°{0,T;V) of the gradient 
system (GS). For this solution, the energy inequality 

(2.5) / \\u\\ 2 H + £ (u(t)) < £ (tio) + / , u) H , teiO.T}, 

Jo Jo 

1 Observe that assuming the compactness of the embedding V > H, as done in |CF10| , is not necessary in our 
case, since the metric on H Y is constant. 
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holds. 

Proof of Theorem \2.7\ The last assertion of Theorem 12.71 follows immediately from the energy 
inequality (|2.5p . because u(t) — u"(t) for all t > 0. It remains to prove that £' maps bounded sets 
into bounded sets and that £ is ^-elliptic. 

To begin with, let / G H Y with < K. By ([22]) we know that ||/|{ ,i}llff xff < CK, 

where C is the same constant as in (12.3[) . Let m be an upper bound for || 0(-) \\hxH on the bounded 
set {y eY : \\h\\hxh < CK}. It follows that 

= sup I f\f'(x) , 6'{x)) H dx + (4>(f\{ ,i}) , 0|{o,i})ffxff| 
\W\ h i,=i Jo 

< ll/llff^ll^llff^ + ||<K/|{0,l})||tfxff||%0,l}IUxH 

< K + mC. 

Finally, in order to show X2-ellipticity we use Lemma \2. 61 Since $ is convex and continuous 
we have 

£x{f) = *(/|{o,i>) > -fc|l/|{o,i}||lxH - P > ~kC 2 \\f\\ 2 CY -(3, fe H Y , 
for some k,/3 £ M. Moreover, 

WfWcy <4||/||^||/|| 2 , feH Y , and H Y ^ C(0; 1; H) ^ X 2 . 
Thus, Lemma [2TB1 yields the claim. □ 

Remark 2.9. The system in Example \2.1\ becomes a gradient system for the choice p\ — P2, 
qi = <72 and X2 — \\q\jp\; furthermore, if \\ < 0, it fulfils the requirements of Theorem \2.1\ and 
possesses a global solution. For A > ; on the contrary, the system undergoes a blow-up in finite 
time as was shown in [ZZ03 . 

Remark 2.10. Observe that if u e H 1 ^, T; X 2 )ni°°(0, T; H Y ), then its C([0,T]; X 2 ) -representative 
satisfies 

(2.6) \\u(t)\\ H i<K 

for some constant K and all (and not only a.e.) t € [0, T]. In fact, consider a Lebesgue-null set 
N on whose complement (|2.6[) holds. Let to £ N and (t n )n£N C (0,T)\N converging to to. Then 
\\u(t n )\\jji < K for all n 6 N and by reflexivity of H Y we deduce that (up to taking a subsequence) 
(w(t n ))neN converges weakly to some v in H Y . It also converges strongly to u(to) in L 2 (0, l;H), 
hence by uniqueness of the limit u(to) = V € H Y . 

3. Invariance properties 

In this section we discuss the issue of invariance properties of solutions to a network gradient 
system. The first results of this kind for solutions of parabolic problems associated with quadratic 
forms go back to Beurling and Deny, cf. [BD59a ( BD59b . A similar criterion in the nonlinear case 
appear already in Brezis' monogr aph (cf. |Bre731 Thm. IV.4.5]), but only recently were Barthelemy, 
Cipriani and Grillo able to prove a more efficient invariance result exactly in the spirit of Beurling 
and Deny. Their main result, which we report for the sake of self-containedness, is the following 
(cf. |Bar961 Thm. 1.1] and |CG031 Thm. 3.4]). 

Theorem 3.1. Let € be a lower semicontinuous, convex functional on a Hilbert space X with val- 
ues in (—00,00], and (T(t)) t >o be the corresponding strongly continuous, nonexpansive semigroup 
on X, generated by the subdifferential £><£. Then (T(t)) t >o leaves invariant a closed and convex 
set C C X if and only if 

<S(P c {x)) < £{x) forallxeX, 
where Pc denotes the orthogonal projection of X onto C. 
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Remark 3.2. Instead of <£(Pc(x)) < £(x) for all x G X , one could of course equivalently require 

Pc{x) G D(<B) and <£(P c (x)) < <£(x) for all x £ £>(<£), 
where D{£) :— {x G X : <£(x) < 00} is the effective domain of €. 

Corollary 3.3. In case the effective domain D{f£) =: D is a Hilbert space in its own right and 
the restriction <&d '■ D — > K is differentiable with derivative <£! D : D — > D' , the condition 

P C (y)€D and <£{P c {y))<<£{y) for ally ED 

is equivalent to 

P C {y)eD and <B' D (P c (y))(y-Pc(y))>0 forallyeD. 

Proof. In both directions the proof is an immediate consequence of 

< Z>(P c (y))(y - P c{ y)) = hm ^(l/) + A( g - P c (y))) - g (P c (y)) < Agfa,) - A g (P c ( y) ) 

a->o A A 

and the fact that P C (Pc(v) + X(y - Pc(y))) = Pc{y)- □ 

Brezis, Barthelemy and Cipriani-Grillo formulate their results in the context of the theory 
of nonlinear semigroups. In fact, under our standing assumptions it is possible to extend the 
Lyapunov function £ defined in (|2.2[) by +00 to the whole X%. Then the extension of £ is lower 
semicontinuous and convex and its subdifferential is single-valued and agrees with the gradient 
V£ (we will denote by £ this extension, too). It is well-known that V£ generates a semigroup of 
Lipschitz-continuous mappings (T(i))t>o on X-z- If € Hy anc ^ V' = 0, then u(t) = (T(t))uo)t>o 
is the solution to 



(GS 



X2 J 



ii(t) +V£(u(t)) = 0, t > 0, 
u(0) = u , 



yielded by Theorem 12.71 Standard references for nonlinear semigroup theory include |Bre73] 
and |Miy92| . In particular, invariance of a closed convex set C under the nonlinear semigroup can 
be rephrased by saying that 

if it e C C X 2 , then u(t) E C for all t > 0. 

Then the following holds. 

Theorem 3.4. Let H be a Hilbert lattice. Under the assumptions of Theorem \2.7\ take a closed 
order interval A C H and consider the order interval 

C := {/ e X 2 : f{x) e A for a.e. x € (0, 1)} C X 2 . 

Denote by Pq the orthogonal projection of X 2 onto C and by Pa the orthogonal projection of H x H 
onto A x A. Consider the following assertions. 

(a) For all initial data u(0) £ C D H Y the solution u to (GSx 2 ) satisfies u(t) € C for all t > 0. 

(b) P c f G H Y and £(P c f) < £{f) for all f G H\ . 

(c) Pcf G H Y and £'(P c f)(f - P C f) > for all felly. 

(d) For all initial data £(0) G (A x A) l~l Y the solution to 



(GSff) 



m + mt)) = 0, t>o, 
m = 6, 



satisfies e Ax A for all t > 0. 

(e) PaC G Y and $(P A C) < *(0 for all £ G Y. 

(f) PaC G Y and $'(PaO(^aO(^ - PaO >0forall£eY. 

Then (a) (6) (c) <= (d) & (e) (/). Moreover, (c) =*> (d) if £ {P c f) = £ a (f) or 
£' (Pcf)(f - Pcf) = for all f G H Y . 
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Proof. The equivalence of (a), (b) and (c) follows from Theorem 13. II Same for the equivalence of 

(d) , (e) and (/), once we consider (GSh) as a (sub)gradient system on the target space H x H . 
Now, observe that 

S(Pcf) = So(Pcf) + Si(Pef), f G H Y . 
We use the fact that, for orthogonal projections Pq onto order intervals of X 2 , Pef G H 1 ^, 1; H) 
if / G H l (0, 1;H) and £q(P c u) < £ (u) (cf. |MuglO| ). Accordingly, a sufficient condition for (6) 
is given by 

PaI\{ ,i}^Y and $(P A /|{o,i}) < *(/|{o,i}) for all / € #f . 
Due to surjectivity of the trace operator 7J 1 (0,l;if) H x H, this condition is equivalent to 

(e) . ' □ 

Example 3.5. We take as order interval on H the positive cone A = H + and accordingly on 
X 2 the positive cone C = X 2 . Assuming Ph+S, = £ + G Y for all £, E Y , Theorem \3.4\ asserts 
that the semigroup (T(t)) t >o associated with (GSx 2 ) is positivity preserving, i.e. T(t)(w ) G X^ 
whenever uq G X% , if and only if $ satisfies < I > (£ + ) < for all £ G Y. The only if part is 

due to £q(/ + )(/ — / + ) = f or a M f £ Hy> ^ s can ^ e seen by identifying the Hilbert lattice H 
with a Lebesgue space L 2 (X) for some finite measure space X , see e.g. [MN91, Cor. 2.7.5], via an 
isometric lattice isomorphism. 

In the nonlinear case, we provide a sufficient condition for well-posedness of a parabolic problems 
in all L p -spaces by applying a Riesz-Thorin-like interpolation theorem. To this aim, it is sufficient 
to show that for any two solutions u,v of (AV) (with initial data uo,vq) the inequality 

\\u(t) - «(t)||oo < ||«o - wollooi *>0, 

holds. This is equivalent to showing that the closed convex sets 

{(f,g)eX 2 xX 2 :\\f-g\\ 00 <a} 

are invariant under the solution to the gradient system associated with the functional 

^(/,s) = #o(/,s) + ^i(/,s) = £(/)+ £(<?), (f,g)eH Y xH Y , 

where 

Mf,9) = £o(f) + £0(9) and #!(/, g) = £ x (f) + £ x {g) = #(/| {0 ,i}) + *(ff|{o,i}). 

More generally, we are led to considering closed convex subsets 

:= {(/,<?) G X 2 x X 2 : (/ - g)(x) G A for a.c. x G (0, 1)} 

of X 2 x X 2 , where A is some closed order interval of H. Similarly, we consider the set 

C A := {(??, 0) G (H x H) x (H x F) : (17 - 6») G A x A}. 

We denote the orthogonal projection of A 2 x A 2 (resp. H x H) onto % (resp. Ca) by P<^ A (resp. 
Pca)- Then the following holds. 

Theorem 3.6. Let H be a Hilbert lattice. Under the assumptions of Theorem \2.7\ take a closed 
order interval A C H . Consider the following assertions. 

(a) For any two initial data (u(0),v(0)) G (H Y x Hy) f~l the corresponding solutions u,v to 
(GSx 2 ) satisfy (u(t), v(t)) G /or aZZ t > 0. 

(b) iV A (/,ff) eH Y xH Y and S{P<^ A (f,g)) < £(f)+£(g) for all f,g€H*. 

(c) Br A (f,g) £H Y xH Y and #'(P % (/, <?))((/, 5) - iV A (/,fl)) > /or aZZ /,<? G 

(cZ,) for anj/ two initial data (£(0), C(0)) G (Y xYJf] Ca the corresponding solutions £, £ to (GSh) 
safe/y £(*)) G C A /or a« i > 0. 

re; Pc A (e, o^xy^ *(Pc A & c)) < *co + *(o v e, c e y. 

(W Pc A (e, C) G y x Y and ^(Po A (f, C))((6 " Pc A (£, 0) > /or aZZ CCeY. 



Of course, such a set agrees with C a setting A := [—a, a]. This makes sense, once we identify the Hilbert 
lattice H with a Lebesgue space L 2 (X) for some finite measure space X, see e.g. IMN91I Cor. 2.7.5], via an isometric 
lattice isomorphism. 
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Then (a) & (b) & (c) <= (d) & (e) & (/). Moreover, (c) (d) if 4>(P VA (f,g)) = #o(/,<?) or 
«SS(iV A (/,ff))((/,fl) - JV A (/,«?)) = /or a /Z / )S 6 4. 

Example 3.7. TTie semigroup (T(t)) t >o associated with (GSx 2 ) * s called order preserving £/ 
T(t)(uo) < T(t)(vo) for all t > and aZZ pairs (tto, «o) G H Y X ify wii/i uq < «o- 
Taking as closed, convex set 

= {(f,g) eX 2 xX 2 :(f- g) £ X+} 
and since (c/. |CG03| Lemma 3.3],) £/ie projection onto < 4>g+ is given by 



f H +{f,9) = I /+ 2 ' ff 2 J 



we can apply Theorem \3.b\ This yields for nonnegative £, cf. [CG03, Theorem 3.8], i/ie sufficient 
and necessary condition 

(3.1) £(f/\g)+£(fVg) <£(f)+£(g) for all f,g € H Y , 

for the semigroup (T(t)) t >o associated with (GSx 2 ) to be order preserving. In our setting it is 
easy to see that actually 

£ (fAg) + £ (fVg) = £ (f) + £ (g) 

whence (|3.ip is eguivalent to 

$(x Ay) + $(x V y) < + for all x, y G Y. 

Example 3.8. In particular, applying the formula for the orthogonal projection onto Ca for 
A := [—a, q]h obtained in [CG03, Lemma 3.3] we obtain the following: 

Let H be a Hilbert lattice. Assume $ and therefore £ to take values in [0, oo), hence V£ to be 
an accretive operator. If 

, / (x — y + a) + (x — y — a)~\ . . , , . / (x — y + a) + (x — y — a)~ 

$ \ x + - — — - - - — -n — - - $ (^) < <%) - $ \y - - — — - + - — k — - 



for all x,y G Y and all a > 0, then the nonlinear semigroup associated with £ is contractive 
with respect to the norms of both X 2 and L°°(0, 1; H). Then by Browder's nonlinear interpolation 
theorem, cf. |GG03[ Theorem 3.6], one concludes that the nonlinear semigroup extends to all spaces 
LP(0,1;H), Kp< oo. 

Example 3.9. Another application of the above invariance criterion shows that if $ is additive 
(resp. homogeneous), then so is the solution operator to (AV): this follows by discussing invariance 
of the closed convex sets 

{(f,g,h):X 2 xX 2 xX 2 :f + g = h} 

( resp. 

{(/,<?) :X 2 xX 2 :af = g} for all aeR) 
under the solutions to the product gradient systems associated with £ . 

4. Dynamic boundary conditions 

We conclude this note by presenting the natural extension of our setting to the framework of 
diffusion equations with dynamic, nonlinear boundary conditions. 

Beside [Bel91, Bcl93, AN93 , we also mention the more recent article |FGGR06j . where a partial 
differential inclusion featuring a linear heat equation and a nonlinear dynamic boundary condition 
is treated. To be more precise, we consider again a closed subspace Y of H x H and discuss the 
Cauchy problem 
(AVD) 

mu(t,x) = £ r u(t,x) + Mt,x), te[o,T], i £(0,1), 

«(%o,i} e Y, te[o,T\, 
&«(%o,i> = -^(^ + 0((«(*)Ko,i}))+*2(*)> te[o,T], * 2 :[0,T]^F, 

u(0) = M , 
W (0)|{0,1} = v . 
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Observe that if Y = Range I for a certain incidence matrix / (see Example 12.21) , ( AVD) is a 
nonlinear generalisation of the setting considered in [MR07] . 

We consider basically the same functional £ introduced in Section [2] but on a different Hilbert 
space. In fact, we introduce the space 



U\ := jf := Q G H\ x H x H : / ]{0jl} - 3 | 



(GS^ 2 ) 



and consider the functional £ : H Y — > R- The reference space is now the Hilbert product space 
X 2 :=L 2 (0,1;H) x Y. 

In this framework the diffusion equation (AVD) corresponds to the gradient system 

u(t) + V£ (u(t)) = #(t), t€[0,T], 
u(0) - u , 

with respect to the Hilbert spaces T-L Y an d ^2 := V2 x Y, where 

"o:=("°) and == , == lM*> *€[0,T]. 

To see this, we assume first f := ( , | G D(V£) and choose a := ( ^ | G with 

UH0.1}/ Wo.i}/ 

5|{o,i} — 0- Thus, we obtain from 

£'(f)fl = / (f(x),9'(x)) H dx + {ct>(f mi} ),0) HxH ± / (&,fl)i, + (6 1 0)j IX fi 
■/o Jo 

that / € H 2 (0,1;H) and /" = — £ and hence for arbitrary g := I J G "H^, after having 

\5|{o,i}/ 

integrated by parts, 

( 1 ff|{0,l}).ffxff + (<H/|{0,1}) 1 9 , |{Q,l})fl"xJ? = (6 , 5|{0,l})ffxH- 

We infer £i=iv(§£ + tf((/|{o,i})). 

For / G H 2 (0, 1; n ffi one can define £0 and £x as above in order to see that f := f , J G 

\/|{o,i}/ 

D(V£). 

Theorem 4.1. Lei : Y — >• Y" &e a function that satisfies 4> = V$ /or some convex <I> G C 1 (Y", R) 
and i/iai maps bounded sets into bounded sets. Then for all initial data u G T-L Y there exists a 
unique solution u G H 1 (0,T;X 2 ) n L°°(0, T; iT|) to (AVD). If * = 0, tten tfie solution also 
belongs to L 2 (0, T; ff 2 (0, 1; i? )) for every T > and even to L 2 (0, 00; H 2 (0, 1; ff)) «/ in addition 
$ and therewith £ is bounded from below. 

Proof. The proof is again based on the application of Theorem 12.81 to the Hilbert spaces H Y 
and X 2 and to the functional £ : H Y K< One sees that £ satisfies (|2.4p . It follows from our 
assumptions that £' maps bounded sets of H Y into bounded sets of (H Y )'. Accordingly, for all 
initial data u G U Y there exists a unique solution u G H 1 (0,T;X 2 ) n L°°(0,T;^) to (AVD). 
Taking the first coordinate of u yields the claim. □ 
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